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Abstract 

quantum pseudo-spin model with random spin sizes is introduced to study the effects of charging-energy disorder on the super- 
^ conducting transition in granular superconducting materials. Charging-energy effects result from the small electrical capacitance 
| ; of the grains when the Coulomb charging energy is comparable to the Josephson coupling energy. In the pseudo-spin model, ran- 
5^ domness in the spin size is argued to arise from the inhomogeneous grain-size distribution. For a particular bimodal spin-size dis- 
tribution, the model describes percolating granular superconductors. A mean-field theory is developed to obtain the phase diagram 
I as a function of temperature, average charging energy and disorder. 
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k*" 1. Introduction 

o : 

— , . Quenched disorder is obiquitous in condensed-matter 
CO systems, sometimes determining entirely new physical 
j» — properties and phenomena. For the theoretical explana- 
C^) tion of many experimental findings, models of classical 
ON and quantum spin systems are employed and - typically - 
disorder is introduced through a spatial variation of the ex- 
change couplings, of the direction of the axial anisotropies, 
and of the sign of the couplings. Theoretical and exper- 
fy^j imental studies of model and real physical systems with 
C3 these kinds of disorder have lead to much new physics. 
Much discussed have been the glassy phases in magnetic 
|1I2I3| , and superconducting materials [1] . 

In superconductors, specially with the advent of high-T c 
ceramic superconductors, the role of disorder has become 
central in the discussion of the physical properties of real 
materials [3|4j . Many materials are naturally microstruc- 
tured and granularity characterises the mesoscopic struc- 
tures of most systems, leading to a phase diagram often dis- 
playing a superconductor-insulator transition at zero tem- 
perature due to the charging energy of the grains. As first 
pointed out by Abeles [5] , when the grain charging energy 
arising from the charge Q and capacitance C of the grain, 
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E c = Q 2 /2C ~ (2e) 2 /d {d is the grain diameter and e 
the electronic charge) is larger than the Josephson-coupling 
energy E a between nearest neighbors grains, phase coher- 
ence is destroyed due to zero-point quantum fluctuations 
of the local phase of the superconducting order parameter 
[5] . Granularity can also be realized in a controlled manner 
in artificially fabricated Josephson-junction arrays of cou- 
pled superconducting grains [7] , with a space dimensional- 
ity d less than 3. Such granular systems can be theoretically 
modelled by pseudo-spin systems, the 'spins' representing 
with their states the few relevant quantum charge states of 
the grains at relatively low temperatures. In such models 
the pseudo-spin size is the same for all grains corresponding 
to the assumption that the grain-size distribution is very 
narrow, or alternatively, of negligible charging-energy dis- 
order. A well known model of this type is the pseudo-spin- 
one model introduced by de Gennes and studied in differ- 
ent works |3|8|9|10j where only charge states —1,0,1 are al- 
lowed, corresponding to S = 1. However, realistic systems 
may contain different kinds of disorder, such as a spatial 
distribution of Josephson couplings between grains or/and 
a distribution of grain sizes, which leads to disorder in the 
grain electrical capacitances and charging energies. 

Studies of the effects of disorder in the electrical capac- 
itance or charging energy of the grains have appeared re- 
cently [11112113] , Within a mean- field approximation [IT] , 
charging energy disorder widens the extent of the super- 



conducting phase at the expense of the insulating one. This 
behavior is consistent with results for the superconductor- 
insulator transition from quantum Monte Carlo simula- 
tions in two-dimensional models of Josephson-junction ar- 
rays [13] with disorder in the diagonal capacitance matrix. 
Earlier calculations for the related boson Hubbard model 
with disorder in the onsite Coulomb repulsion [TJ] are also 
consistent with a decrease in the phase coherence thresh- 
old and moreover suggest a different universality class from 
the non-disordered case [15] . 

In this work, we consider the effects of charging energy 
disorder in granular superconducting materials within the 
pseudo-spin approach. By considering a spatial distribu- 
tion of the grain sizes, which leads to local charging en- 
ergy disorder, a quantum pseudo-spin model with random 
on-site spin sizes can be constructed. A mean-field theory 
is developed to obtain the phase diagram as a function of 
temperature, average charging energy and disorder. 

Spin-size disorder models have seldom been considered in 
the literature. In the context of classical spin models, spin- 
size disorder can be readily turned into exchange-coupling 
disorder; the physics of quantum systems with spin-size 
disorder, however, appears to have not been investigated 
in such depths. It has been considered mostly within the 
one-dimensional systems: dilution of a quantum spin-i 2- 
ladder was studied by Sigrist and Furusaki lj whilst the 
general problem of a quantum spin chain with random S as 
well as random J was considered by Westerberg et. al [5] 
within a real-space renormalization-group method showing 
that these systems belong to a different universality class 
of disordered spin systems. 

2. Pseudo-spin model with random spin sizes 

The Hamiltonian for a set of superconducting grains cou- 
pled by the Josephson energy can be written as the sum of 
the Coulomb charging energy and the Josephson-coupling 
energy |3I6I8| 

H 9s = ~ E G V l QiQi - E E v cos & ~ *i)> 

i,j <ij> 

where Qi is the net charge on the superconducting grain at 
site i and 9i is the phase of the local superconducting order 
parameter fy. is the electrical capacitance matrix and 
Eij is the Josephson coupling between nearest-neighbor 
grains. The charge Qi in each grain can be expressed in 
terms of the excess number of Cooper pairs m as Qi — lern. 
Considering a diagonal capacitance matrix Cy = CiSij and 
uniform Josephson-coupling energies Eij = E Q , leads to the 
self-charging model 

H gs = 2 J2 u i n i - E costfi - 6j), (2) 

i <ij> 

where Ui = e 2 /Ci are the charging energies of the grains. 
The number operator n, is conjugate to the phase 9i, sat- 
isfying the commutation relation 



[m,9j} = -iSij, (3) 

and can be written as — —id/dOi, having integer 
eigenvalues 0,±1,±2.... The model of Eq. ((2|) can also 
be regarded as a boson Hubbard model |14|15j where the 
charging energy represents the onsite Coulomb repulsion of 
bosons and the Josephson coupling represents the hopping 
term. 

When the charging energies are uniform Ui — U , a 
pseudo-spin one model for the Hamiltonian of Eq. ([2|) can 
be constructed by truncating the basis vectors of the num- 
ber operators |n$ > to |0 > and | ± 1 >, corresponding to 
the lowest charging energy states. Identifying the charging 
states, | n$ >, as the eigenstates of S z for spin S = 1, writ- 
ting the second term of Eq. (J2j) in terms of e ±l6i and making 
the correspondence 

e i9 < -> S+/V2, e~ ie i - SJ/V2, n z -» S?, (4) 

one obtains the de Gennes mapping to a S = 1 pseudo-spin 
model with single-ion spin-anisotropy |3l8j 

Hs=i = 2U a Y,(SD 2 ^ E(^ +5 7 + S t S ^)> ( 5 ) 

with S+ = Sf + iSf. The factor l/y/2 in front of S± in 
Eq. (j4|) takes into account the length of the pseudo spin 
y S(S + 1). In terms of this spin model, the superconduc- 
tor to insulator transition in the original granular supercon- 
ductor corresponds to a zero temperature transition where 
the ferromagnetic order of the spins in the xy plane is de- 
stroyed by quantum fluctuations when the S z component is 
confined to zero for an increasing ratio of U /E . Although 
this mapping is not exact [18j . the T — critical behav- 
ior observed in the original phase model of Eq. without 
disorder and the spin model is the same [14119] as found 
in numerical calculations in one dimension. For higher di- 
mensions, results from the same mean-field approximation 
applied to both models also agree [8] . 

To consider the main effects of disorder in the charging 
energies Ui, we generalize the above approximate mapping 
to a pseudo-spin model with randomness in the spin-size 
values Si. Since for a given charging energy (or tempera- 
ture fluctuation), low values of Ui correspond to charging 
states with higher rij, it seems reasonable to use a trunca- 
tion scheme which identifies the charging states |rij > and 
\n,j > with the S^-eigenstates of the spins of different sizes 
Si and Sj , such that the corresponding maximum charging 
energies 2UiSf and 2UjSj of the grains are both compara- 
ble to the same truncation energy (here Si = \J Si(Si + 1) 
is the length of the spin S.^). This (approximate) mapping 
leads to the effective spin Hamiltonian 

H = D E i(^) 2 + J E tV^7 + sfsn, (6) 

with randomness in the spin-size values Si. Here, J = 
—E a /2, D = 2US 2 where U is the average value of Ui and 



S the average value of Si. The spin sizes Si are restricted 
to take only integer values 0, 1, 2, 3, ... in this mapping. In 
Eq. ([6]) , randomness in the spin-size also leads to random- 
ness in the local single-ion anisotropy parameter, but these 
random variables are only correlated at the same site. 
A typical spin-size distribution we shall be considering is 



P(Si 



: E 

L=±l 



5{Si - S - L) + (1 - 2x)S(S i - S ), (7) 



where three spin-size values S , 5 G + 1 and S — l are mixed 
with average value S and a concentration x, which is a 
measure of the disorder. It is also convenient to rewrite the 
spin model of Eq. ([6]) in terms of normalized spins Si = 

Si /Si 



h = dJ2(s[) 2 + -/X :S V' S . + § t § n- 

* («■> 



(8) 



In this latter form, a bimodal distribution of spin values 
can be considered: 



P{Si) = x6(§i -&) + (!_ x )6(Si - 5 2 ), 



(9) 



which in the special case S\ = corresponds to the dilu- 
tion of a spin-S2 system or to percolating granular super- 
conductors [TB] in the presence of charging effects. 

3. Mean-Field Theory 

For high enough space dimensions d, reasonable results 
can be obtained from simple mean-field approximations 
(MFA). The MFA replaces all variables around a given 
small cluster with their average value, or order parame- 
ter, and then use the resulting approximate Hamiltonian to 
evaluate the order parameter itself. Since the spin sizes are 
random, one must carry out the averaging procedure with 
respect to P(Si) as well and we denote this by [...] ace . 

The simplest cluster is a single spin on a site. With M = 
[< Si >] WV e the mean field, this yields a self-consistent 
equation for the order parameter M: 



M = [(Si 



1 



-'MFA 



-Tr s ,Sie 



-0H M1 



(10) 



where Hmfa is the mean-field Hamiltonian and (3 = 
1/ksT. For the anisotropic Hamiltonian of Eq. ©, we 
assume that the ferromagnetic ordering takes place in the 
xy plane of the anisotropy, say [< Sf >] a ve — M. Then, in 
the MFA we replace [3] 



(S?)S? + Sf(S?)-(S?)(Sf), 



(11) 



OX ni I aV qV 
D i D j "T D i D j 

by virtue of the fact that < Sf >— 0. Near the critical 
curve, the MFA Hamiltonian then becomes: 



Hmfa 



\S x ) + 0(M 2 ), 



(12) 



where A = 2z\J\M/D is a dimensionless expansion param- 
eter. 



Since we are interested in the phase transition line only, 
we can assume M = M/S small and use a first-order per- 
turbation expansion to solve Eq. (|10p for arbitrary spin S. 
In a now single-site problem, we reabsorb a factor S~ 2 into 
D and J. The unperturbed states |m > are eigenstates of 
S z and thus also of the unperturbed Hamiltonian H Q = 
D{S Z ) 2 , with to = -5, -S + 1, S-1,S. We have 



(S x ) = 



E, 



-/3E„ 



(lp m \S X \lpr, 



Er 



-0E„ 



where to first-order E„ 
and 



\lprn) = | TO) + 



m' ,m 



Dm 



XDS X \ 



(13) 

(ml -XDS x \m) = Dm 2 



E° - E° , 



(14) 



Here, the prime on means ml ^ to. These can then be 
inserted in Eq. (fT3"|) to get, keeping only first order terms 
in A: 



(S x ) 



Zo 



where Z\ 



o 



£ 



E 



_ 0Dm , (m\S*\m')(m'\S x \m) 



e -l3Dm 2 _ Uging 



A. 



(15) 



(m\S x \m') = + A s m 5 m , tm+l ) 



Ai = ^(S-m)(S + m+l), 



wc get 

(S*) = 



(A-m-l) 2 



2Z (to + l) 2 - TO 2 ^ (to - l) 2 - TO 2 



(16) 



,-(3Dm 2 



which, using Eq. (jlOp . gives implicitly the critical temper- 
ature in the MFA 



D 



= [ 



1 



2z\J\ E m rDm!!/(S(5+1)terc ) 



E 

m=-S 



5(5 + 1) 



1 — 4to 2 



,-_Dm 2 /(S(S+l)fcsT c )] 



,(18) 



One can check that for 5 = 1 and in absence of disor- 
der the MFA transition line from Eq. (fT8|) agrees with the 
result obtained for the de Gennes model of Eq. ((5]) in the 
same approximation [3] . Fig. 1 shows the phase boundaries 
between the superconducting and normal phases obtained 
in absence of disorder for different values of uniform spin 
sizes 5 = 1,5 = 2 and 5 = 5. The critical value E /U for 
the superconductor-insulator transition at T = remains 
unchanged and only at much higher temperatures there is 
a deviation in the transition lines when considering larger 
spin approximations. This is in agreement with MF calcu- 
lations using truncated basis vectors in the phase model of 
Eq. 10 for increasing number of states [8] . 
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Fig. 1. Phase boundary without disorder obtained from the pseu- 
do-spin model with different uniform spin values Si = 1,2 and Si = 5. 

Our present calculations generalize the MFA to arbitrary 
spin sizes and allow to obtain the transition line in the case 
of a random spin size distribution by simply averaging with 
respect to P(S) in Eq. (fTS]) . The consequences of spin-size 
disorder depend very much now on the chosen distribution, 
but it is expected that in general within a simple MF ap- 
proximation, the phase transition of the homogeneous sys- 
tem is preserved with at most a modification of the phase 
boundary line. Fig. 2a shows the transition lines for the 
spin-size distribution of Eq. (|7|) with S = 2 and different 
values of the disorder parameter x. For increasing disor- 
der the critical value for phase coherence E /U decreases, 
increasing the extent of the superconducting phase at low 
enough temperatures. This behavior is in qualitative agree- 
ment with MF calculations [11] in the phase model of Eq. 

On the other hand, a different behavior is expected when 
there is dilution of superconducting grains. Fig. 2b shows 
the transition lines for the spin-size distribution of Eq. ([9]) 
with Si = 0, S2 = 2, corresponding to a dilution of grains 
for different concentrations x. In this case, disorder de- 
creases the extent of the superconducting phase. For large 
values of x, corresponding to the system below the perco- 
lation threshold, the superconducing phase should disap- 
pear. However, an improved MF approximation is required 
to describe this behavior. One possible approach is the MF 
renormalizaton-group method [17j used previously to study 
the de Gennes model of Eq. §5§ in absence of disorder [9] . 

4. Conclusions 

We have introduced a quantum pseudo-spin model with 
random spin sizes to model the effects of charging-energy 
disorder in granular superconducting materials. Random- 
ness in the spin size is argued to arise from the inhomoge- 
neous grain-size distribution. For a particular bimodal spin- 
size distribution, the model describes percolating granu- 
lar superconductors. A mean-field theory has been devel- 
oped to obtain the phase diagram as a function of tempera- 
ture, average charging energy and disorder. The results are 
qualitatively consistent with previous mean-field calcula- 
tions in the phase-number representation. The pseudo-spin 
model should provide a useful framework to study the crit- 
ical behavior and universality classes in presence of strong 
charging-energy disorder. 
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Fig. 2. Phase boundary with disorder obtained from the pseudo-spin 
model with random spin sizes: a) with charging energy disorder, 
corresponding to disordered spins values giving by the probability 
distribution of Eq. (0, with average spin value So = 2 and different 
x ; b) with dilution of grains, corresponding to the disordered spin 
values with the probability distribution of Eq. J9j, with S± = and 
S2 = 2 and different x 
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